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Abstract. 

We calculate the Euclidean action of a pair of Z 2 monopoles (instantons), as a function 
of their spatial separation, in D=2+l SU(2) lattice gauge theory. We do so both above 
and below the deconfining transition at T = T c . At high T, and at large separation, we 
find that the monopole 'interaction' grows linearly with distance: the flux between the 
monopoles forms a flux tube (exactly like a finite portion of a Z2 domain wall) so that 
the monopoles are linearly confined. At short distances the interaction is well described 
by a Coulomb interaction with, at most, a very small screening mass, possibly equal 
to the Debye electric screening mass. At low T the interaction can be described by 
a simple screened Coulomb (i.e. Yukawa) interaction with a screening mass that can 
be interpreted as the mass of a 'constituent gluon'. None of this is unexpected, but it 
helps to resolve some apparent controversies in the recent literature. 
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1 Introduction 



In a recent paper [[!]] the free energy of a separated pair of Z 2 monopoles has been 
calculated in D=3+l SU(2) lattice gauge theory. This problem is difficult because it 
involves the numerical evaluation of the ratio of two different partition functions, so that 
the calculations in |j] are only accurate at small monopole separation. The conclusion 
is that the interaction potential has a simple Coulomb/Yukawa form both at high and 
at low T; and that this remains so, at high T, even if one considers monopoles whose 
world-lines are spacelike rather than timelike. While the other results are consistent 
with expectations, this last is not. What one expects is that the Z 2 flux between the 
monopoles will form a 'flux tube' so that the interaction energy grows linearly at large 
separation. Motivated by this unexpected result we have performed similar calculations 
using a different method that is less direct but much more accurate. This method has 
been used recently [0, [J in closely related investigations of vortices || and, some time 
ago, of high-T Z 2 domain walls ||. Our calculations are carried out in 2+1 Euclidean 
dimensions where, apart from the fact that the monopoles are instanton-like rather than 
soliton-like, not only can the same questions be posed as in 3+1 Euclidean dimensions 
but they can be answered much more accurately. In addition the close relationship 
between Z 2 monopoles, vortices and domain walls is particularly transparent in 2+1 
dimensions. (Although it holds in D=3+l equally well.) As we shall emphasise, after a 
suitable relabelling of the space-time axes, one can see |2], |J that the high-T Z 2 domain 
wall is a periodic vortex line (as created by a periodic 't Hooft disorder loop ||) since 
this phase is a Higgs phase for vortices. And in this phase the Z 2 flux between a pair 
of well-separated monopoles passes through a flux tube that is nothing but a portion 
of such a domain wall that ends on these monopoles; so that such Z 2 monopoles are 
'linearly confined'. We shall indeed confirm that the interaction is linear at high T, 
as expected. The contrary conclusion in jj]] can presumably be attributed to the lack 
of accurate calculations at larger inter-monopole separation. We have not pursued the 
same calculations in D=3+l since, as our D=2+l calculations were completed, a paper 
appeared || addressing similar questions in D=3+l with conclusions similar to ours 
(where they overlap). 

In the next section we describe Z 2 monopoles and explain how we expect them to be- 
have. We then describe our Monte Carlo calculations and finish with some conclusions. 
Since most of the techniques and physical quantities involved have been described in 
detail in earlier papers |l|, |2|, [5|, we shall aim to be succinct rather than pedagogical in 
this brief report. 

2 Z2 monopoles - expectations 

In 2+1 Euclidean dimensions, the monopoles are localised in space and time and so 
can equally well be thought of as instantons. What we shall do is calculate the extra 
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Euclidean action, AS, associated with the presence of a pair of such monopoles. If we 
were in D=3+l this quantity would be proportional to the derivative of the free energy 
of the monopole pair, and its variation with the monopole distance would tell us how 
the monopole interaction energy varies with distance. In D=2+l we shall, by analogy, 
sometimes refer to the action difference as the monopole 'interaction energy'. 

To calculate the interaction between Z 2 monopoles at sites x and x' (on the dual 
lattice) one introduces a Z 2 Dirac string that starts at x and ends at x' . Such a string 
flips the sign of any plaquette that it threads. (This is so if, as here and in the related 
literature, the SU(2) group elements on the links are represented by fundamental 2x2 
matrices rather than by adjoint 3x3 matrices.) One can introduce such a Dirac string 
by changing j3 — > —(3 on the plaquettes threaded by the string. We shall refer to this 
as a set of 'twisted' plaquettes because, in the case where one obtains x' by translating 
x by a single lattice period, this is equivalent to the imposition of twisted boundary 
conditions. By redefining suitable subsets of link matrices by U' — — U fl], || [5[], it is 
easy to see that the actual location of the string, apart from its end points, is invisible; 
as we would wish for a Dirac string. 

2.1 T > T c 

Let us first consider the extreme case where the line of twisted plaquettes extends right 
across the lattice; say in the x-direction with the twisted plaquettes in the (y, t) plane. 
This produces a Z 2 Dirac string that closes upon itself through the boundary; so it 
is not attached to any monopoles. Such a twist is, as remarked earlier, equivalent to 
imposing twisted boundary conditions. For T > T c it imposes the presence of a Z 2 
domain wall separating the two high-T Z 2 vacua ||. As pointed out in 0, ||, this 
domain wall is in fact a periodic spacelike vortex tube that has been squashed by the 
narrow extent in the Euclidean time direction. Its existence demonstrates that such 
spacelike Z 2 flux is, at high T, confined to a flux tube. So we expect that for T > T c 
the interaction energy between two spacelike separated Z 2 monopoles will grow linearly 
with their separation, r, once this separation is larger than the width of the tube 

AS ~ ar r large. (1) 

(Note that all quantities will be expressed in lattice units from now on.) At small 
separations the flux of one monopole spreads out in the three Euclidean space-time 
directions and then returns to the other monopole with no opportunity to form a flux 
tube. Thus the interaction should be Coulombic cx 1/r. As r is increased we might 
expect a Yukawa-type correction 

AS" cx r small. (2) 

r 

with the scale /1 being plausibly related to the Debye electric screening mass [[J. The 
latter governs the width of the vortex tube/domain wall which, for large r, carries the 
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return flux between the monopoles over most of the distance between them. 

Note that if there were no linear term then the asymmetric geometry of our system 
would mean that, as we increased r to r > 1/T, any 1/r Coulomb term would turn 
into a logarithm oc log(rT). This is also the case if we consider spatially separated Z 2 
monopoles at high T in 3+1 dimensions. We shall ignore this for the purposes of our 
fits because, in practice, the Coulomb or Yukawa contribution becomes very small for 
r > 1/T. 

At very high T we expect the tension of the Z 2 domain wall to be calculable in 
perturbation theory and, as shown in , the perturbative expression is accurate down 
to remarkably small values of T/g 2 . (Recall that in D=2+l g 2 has dimensions of mass.) 
Thus the fact that our spacelike separated monopoles are joined by a section of such a 
domain wall, allows us to estimate the value of a in eqn(|l|): 

3 

where (3 is the usual (inverse) lattice coupling and 

a = a pert (l + 0.25^) (4) 

provides an accurate representation of the Monte Carlo calculated domain wall tension 
for g 2 /T < 1.2. The leading order perturbative value, a pert , can be calculated as a 
function of the temporal lattice discretisation and the resulting values are listed in 
Table 8 of §. 

In addition to the above quantitative expectation for the interaction between spa- 
tially separated monopoles at high T, one might equally well consider the interaction 
of two monopoles with a timelike separation. However, because such a separation is 
limited to very small values, r < 1/2T, there is nothing very interesting to be learned 
and we shall not consider that case any further here. 

2.2 T < T c 

Consider now the low-T confining phase. In this case if there are vortex tubes, they 
condense into the vacuum and thus the interaction between monopoles becomes inde- 
pendent of r at large r. At small r we again expect a screened Coulomb interaction 
oc exp(— /J,r)/r, or a sum of such terms. But now there is no T to set the scale, so 
(the lightest) /i will presumably be related to the mass gap, which is provided by the 
lightest ++ glueball in this theory 0. 

The above discussion translates directly to D=3+l. Now the Z 2 monopoles will 
have world lines and the Dirac string will sweep out a sheet, which requires a sheet of 
twisted plaquettes. If we are at high T and the monopole world lines and separations 
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are spacelike, then the return flux once again goes through flux tubes, with the scale 
set by the electric Debye mass. So at large separations the interaction energy should 
grow linearly and at small distances it will be screened Coulombic. By contrast, at low 
T the vortex tubes condense in the vacuum and there is only a screened Coulombic 
interaction. While the calculations of agree with the above expectations at low T 
they claim that a simple Coulomb interaction, with no linear piece, works at high T 
as well. Our method of using the action difference is, however, much easier and more 
accurate, especially at larger r, than the method employed in |l| and we now turn to 
calculations designed to resolve this puzzle. 

3 Z2 monopoles - calculations 

Our Monte Carlo calculations use the usual plaquette action parameterised by the 
(inverse) coupling (3 where (3 — > A/ag 2 as a — > 0. The low-T physics of this theory is 
discussed in ||, the properties at high T in ||, and the deconfining transition in [[|. 

We perform calculations at two values of the coupling: (3 = 9 and (3 = 13.5. In each 
case we perform one calculation at T > T c and one at T < T c . At (3 = 9 we perform 
calculations on 24 2 4 and 24 2 8 lattices and at (3 = 13.5 on 36 2 6 and 36 2 12 lattices. 
On a L 2 L t lattice of spatial extent L ^> L t we have aT = 1/L t so the above choices 
correspond to higher and lower T respectively. To be more quantitative we recall that 
the critical value of L t that corresponds to T = T c is given, as a function of (3, by || 

^09)^^03-0.37). (5) 

We see that, at both values of /?, the lower value of T corresponds to T ~ 0.7T C while 
the higher value corresponds to T ~ 1.4T C . 

For each of these lattices we calculate the average action with a line of twisted pla- 
quettes of length r which introduces two Z 2 monopoles a distance r apart as discussed 
earlier. For each value of r we perform a Monte Carlo calculation, typically with 10 6 
sweeps, of the average total action, (S(r)). The action here is S = Z) p (l — ^TrU p ) 
where U p is the product of the four SU(2) matrices around the plaquette p, and c p is 
— 1 for the twisted plaquettes and +1 otherwise. We plot in Tables || and @ our results 
for the action. We turn now to a discussion of these results. 

3.1 T > T c 

In Fig.|l| we plot 

AS = (S(r)) - (S(0)) (6) 

for the (3 = 13.5 high-T calculation. We observe an unambiguous confirmation of the 
expected linear rise with r at large r. Note that the point at r = 36 has not been 
misplaced. Its value is given by a periodic flux tube with no monopoles attached; thus 
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it is missing the large (divergent as a — > 0) Coulombic self-energy of the monopoles, 
that otherwise contributes to S{r) (except when r is close to the lattice period). 
We parameterise AS" using a sum of linear and Yukawa terms 

AS = c + c h or. (7) 

r 

This conventional parameterisation is clearly approximate. We do not expect the 
linear piece to be there at all at small values of r. In addition, as remarked earlier, the 
Coulomb 1/r should mutate into a logr dependence at larger r. Moreover we expect 
a string/roughening correction to the linear term at large r, which can be calculated 
in perturbation theory at high T ||. All of these corrections to eqn(^) are small and 
we shall assume, perhaps optimistically, that they are negligible for our purposes. 

We have performed fits with eqn(|7|) to, for example, the range r £ [1, 30]. We exclude 
r = from the fit for obvious reasons and we also exclude values r > 30 since as r 
approaches the lattice period, at some point it becomes energetically favourable for the 
'Coulombic' part of the monopole-monopole flux to join through the periodic boundary 
and for this to be matched by a periodic flux tube whose location is now decoupled 
from the location of the twist. We find that fits as in eqn(^) possess a good \ 2 ~ about 
0.9 per degree of freedom for the best fit shown in Fig.[l|. There are two immediate 
questions: is the monopole string tension as expected and what is the screening mass? 

First, consider the monopole string tension. From the fit in eqn(^) we obtain a value 

a = 0.0300 ±0.0009 (8) 

which is consistent with the value obtained from the r = 36 value: a = 0.0315(12). 
We can compare this with our high-T theoretical expectation in eqn(|3]). Plugging in 
(3 = 13.5, T/g 2 = (3/AL t = 13.5/24, and the value of a from eqn(^) using the L t = 6 
value of a pert as listed in Table 8 of we obtain an expected value a ~ 0.035. If 
instead of employing the heuristic correction term in eqn(§) we use only the one- loop 
value for a then we obtain a ~ 0.024. Thus our fitted value of a is in good agreement 
with our expectations; remarkably so given that the one-loop perturbative expression 
has been derived for much higher T than the values being considered here. 

A similar analysis of the (3 = 9 calculation on a 24 2 4 lattice yields best fits with 
a poorer, albeit not unacceptable, x 2 - We obtain a monopole string tension a = 
0.078(3) which is consistent with the string tension one finds with a periodic twist, 
a = 0.0735(15). This can be compared with what one obtains from high-T perturbation 
theory, eqns([|): a ~ 0.056 without the heuristic 0(g 2 /T) correction in eqn(^), and 
o ~ 0.081 if we include it. Satisfactory agreement, in other words. 

We turn now to the screening mass fi. Its fitted value is: 

fj. = 0.000 ± 0.044 /? = 13.5 (9) 

That is to say, we are quite consistent with a zero mass and so a purely Coulombic 
(unscreened) short distance interaction between the monopoles. Certainly the screening 
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mass is much less than the T = mass gap. At high T there are additional natural mass 
scales, the simplest being aT c ~ 0.12, aT = 1/6, \J g' 2 T ~ 0.22. All of these are probably 
too large to be consistent with our value of \x. There are other more dynamical mass 
scales such as the (electric) Debye screening mass, m D , and the spatial string tension 
a s . We have performed an explicit calculation of Polyakov loop correlations on the 
same 36 2 6 lattice at (3 = 13.5 in order to determine these masses ||, and we find 

m D = 0.103 ±0.005 (3 = 13.5 (10) 

and 

y/as = 0.139 ±0.002 /3 = 13.5. (11) 

The value of y/a s is too large for fi but there is perhaps a possibility, at the 2a level, 
that the Debye screening mass provides this scale. 
A similar analysis of the (3 = 9 calculation yields 

/i = 0.2 ±0.2 ^ = 9.0. (12) 

This is again consistent with the (electric) Debye screening mass that we have calculated 
on the same lattice at the same (3, 

m D = 0.16 ±0.01 (3 = 9.0 (13) 

although, because of the larger errors, it is also consistent with most of the high-T 
scales referred to earlier (but not to the T = mass gap). 

The value of fi is hard to determine numerically because, as we see from Fig.l, even 
at our relatively large value of (3 our resolution of the region in which the 'Coulombic' 
interaction dominates is poor. We have tried fitting the values at small r with a purely 
Yukawa form, and that again produces a value consistent with zero. We have also 
modified the form in eqn(^) by multiplying the linear piece by a factor 1 — exp(— /iV). 
This is quite natural; after all one would expect the linear term to gradually disappear 
once r becomes less than the flux tube width. Such fits give values of \x and \J that 
are non-zero and consistent with the value of in egnflTop. 

3.2 T < T c 

In Fig.|^ we plot the action, AS, of the two Z 2 monopoles as a function of their sepa- 
ration, r, in the low-T confining phase at (3 = 13.5. In contrast to the situation in the 
deconfined high-T phase, we see that there is no linearly confining interaction between 
the monopoles. This is as expected, spacelike Z 2 vortex tubes disorder timelike Wilson 
loops, and, in the confining phase, will condense in the vacuum, thus costing no extra 
action. 
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At small separations the monopoles feel a Coulomb interaction which should be 
screened at larger r in this phase. Thus we fit AS to the Yukawa form: 



AS = c + c . (14) 

r 

We find that an unscreened Coulomb interaction is statistically unacceptable (see the 
fits shown in Fig.|2|); there is clearly a non-zero screening mass 

fi = 0.4±0.1 /? = 13.5. (15) 

Similar low-T calculations on the 24 2 8 lattice at (3 = 9 also reveal a zero string tension 
and a non-zero screening mass: 

/i = 0.43 ±0.09 p = 9.0. (16) 

Since the above masses are in lattice units one would expect the f3 = 13.5 mass to 
be about two-thirds of the (3 = 9 value; which is quite possible within the rather 



large errors. We also observe that taken together the values in eqns([L5],[T6p are not 
really compatible with the T = mass gap f|: mo — 0.495(10) at (3 = 13.5 and 
me — 0.76(1) at (3 = 9. This is reminiscent of what has been found JT0] for the 
screening of the Abelian monopole flux in the maximally Abelian gauge in D=3+l. 
The screening mass there was interpreted as the constituent gluon mass: half the 
average mass of the lightest + and 2 + glueballs. We note that our screening masses, 
in eqns(0,|IBp, are consistent |3J with being interpreted in this way. 

4 Conclusions 

We have shown that in the low-T confining phase the interaction of a pair of Z 2 
monopoles a distance r apart can be described by a simple Yukawa interaction cx 
-exp(— fir) where the effective screening mass \i is consistent with what one might 



expect for the mass of a 'constituent gluon' flQfl. In the high-T deconfined phase the Z 2 



flux between widely separated monopoles is collimated in a flux tube and the interaction 
energy increases linearly with distance, i.e. the monopoles are 'linearly confined'. Such 
a spacelike flux tube, when closed upon itself through a spatial boundary, becomes a 
domain wall separating two high-T Z 2 vacua. Its string tension can be calculated in 
perturbation theory at very high T, since the dimensionless coupling, g 2 /T, becomes 
small there. The string tension we obtain is surprisingly close to the value expected 
from one- loop perturbation theory, given that in our 'high-T' calculations g 2 /T is in 
fact close to unity. At small r, where the flux tube has not yet formed, we find that the 
interaction is consistent with a Coulomb interaction possibly modified by very weak 
screening. This small screening mass is just about consistent with the Debye electric 
screening mass that we have also calculated. All these calculations are for values of the 
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lattice spacing a that are sufficiently small, expressed either in units of g 2 or in units 
of T, that we can expect negligible corrections when taking the continuum limit. 

All our calculations have been for monopoles in 2+1 Euclidean dimensions where 
monopoles are instantons rather than solitons. However we expect things to be very 
much the same in 3+1 Euclidean dimensions. At high T there are, again, perturbatively 
calculable 'domain walls' between differing Z 2 vacua, and these walls are just the regions 
swept out by a periodic Z 2 flux tube as it moves through space-time in a spacelike 
direction. Consider two Z 2 monopoles joined by a Dirac string made invisible through 
the use of appropriately twisted plaquettes. If their spacelike separation is large and 
they propagate in a spacelike direction, the Z 2 flux between them will generate a 
portion of such a Z 2 domain wall, that is bounded by the monopole world lines. If we 
relabel the direction of propagation as the time direction, so that the short Euclidean 
t direction becomes a spatial direction, then we can re-express the above as telling 
us that, within this three dimensional space in which one of the spatial directions is 
short (and equal to 1/T) but the temperature is now zero, the potential energy of two 
Z 2 monopole sources is linearly confining. (The asymmetric spatial geometry will by 
itself enforce confinement, if only logarithmic, at large r.) Moreover the corresponding 
string tension is calculable in perturbation theory. At low T by contrast, we expect 
any such flux tubes to condense in the vacuum because the system is confining, so that 
the monopoles experience a screened Coulomb interaction just as in 2+1 dimensions. 

The aspect of the present calculations that could be readily improved is that of the 
screening masses. This would require calculations at smaller a so as to increase the 
resolution of the region in which the Coulomb/Yukawa potential is visible. If one were 
to consider not just the total action, but the action over various distances around each 
or both of the monopoles, the statistical accuracy of the calculations would be greatly 
increased as well. 

Note added As we were preparing to submit this paper for publication, a revised 
version of [M appeared in which a linear fit to the high-T inter-monopole interaction is 
used. This was shown to be necessary in the recent, more accurate, D=3+l calculations 
of ||, and, as we have shown in this paper, the same is true in D=2+l. 
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Table 1: The total lattice action, S, on 24 2 4 and 24 2 8 lattices at (3 = 9, as a function 
of the separation, r, between the two Z 2 monopoles. 
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Table 2: The total lattice action, S, on 36 2 6 and 36 2 12 lattices at /3 — 13.5, as a 
function of the separation, r, between the two Z 2 monopoles. 
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Figure 1: The action of two Z 2 monopoles a spacelike distance r apart on a 36 2 6 
lattice at (3 = 13.5 for which T > T c . 



12 



3 



in 
< 



o 




^ ' . . »-»-r*"»- 



Yukawa fit 

■ Coulomb fit 







8 



12 



16 



20 



Figure 2: The action of two Z 2 monopoles a spacelike distance r apart on a 36 2 12 
lattice at /3 = 13.5 for which T < T c . 
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